Pressure effects on the superconducting transition in nii-CaAISi 
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We present a combined experimental and theoretical study of the effects of pressure on T c of the 
hexagonal layered superconductors nif-CaAlSi (n = 1, 5, 6), where nH labels the different stacking 
variants that were recently discovered. Experimentally, the pressure dependence of T c has been 
investigated by measuring the magnetic susceptibility of single crystals up to 10 kbar. In contrast to 
previous results on poly crystalline samples, single crystals with different stacking sequences display 
different pressure dependences of T c . l//-CaAlSi shows a decrease of T c with pressure, whereas 
5H and 6-ff-CaAlSi exhibit an increase of T c with pressure. Ab-initio calculations for 1H , 5H and 
6H -CaAISi reveal that an ultrasoft phonon branch associated to out-of-plane vibrations of the 
Al-Si layers softens with pressure, leading to a structural instability at high pressures. For 1H- 
CaAlSi the softening is not sufficient to cause an increase of T c , which is consistent with the present 
experiments, but adverse to previous reports. For 5H and 6H the softening provides the mechanism 
to understand the observed increase of T c with pressure. Calculations for hypothetical 2H and 3H 
stacking variants reveal qualitative and quantitative differences. 

PACS numbers: 74.70.Dd, 74.62. Fj, 74.25.Kc, 74.62.-c, 74.25. Jb 



Superconductivity in hexagonal layered compounds 
has attracted broad interest since the discovery of " high- 
ly superconductivity in MgBz 1 -^ and other structurally 
related compounds such as (Ca,Sr)AlSi f 3 - CaSi2^ and 
very recently, alkaline-earth intercalated graphites^ii 
In all these compounds the light elements are arranged 
in honeycomb layers which are intercalated by alkaline 
earth atoms. Depending on the elements forming the 
honeycomb layer, the electron and phonon states in- 
volved in the superconducting pairing are quite differ- 
ent. In MgB2, holes in the a bands of the B layer cou- 
ple strongly to B bond-strctching phonon modes. In 
CaSi2, which can be considered as the "antibonding ana- 
logue" to MgB 2 , the a* bands are strongly coupled to the 
bond-stretching phonons of the Si layer. In CaC6 and 
CaAISi, the so-called "interlayer" bands are filled and 
experience significant electron-phonon (e-ph) interaction 
to the out-of-plane buckling vibrations of the honeycomb 
layers i 
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CaAISi is of particular interest because it exhibits an 
ultrasoft phonon mode and crystalizes with several stack- 
ing variants ! 13 ! 14 The interplay of these two ingredients 
gives rise to intriguing effects both on the normal and the 
superconducting properties. The presence of an ultra- 
soft phonon mode at ~7 meV was initially evidenced by 
first-principles calculation a 11 ' 12 ! 15 ' 16 : 17 and recently con- 
firmed by neutron-scattering experiments^ The ultra- 
soft phonon modes are believed to induce strong e- 
ph coupling and to lead to an enhanced specific heat 
anomaly at T c as well as to a positive pressure depen- 
dence of T c in CaAISi. These effects have not been ob- 
served in the isoelectronic and isostructural SrAISi which 
does not display any signature of soft phonon modest 




FIG. 1: (color online) Crystal structure of 1H, 5H and 6H- 
CaAlSi, which are characterized by a different sequence of A- 
and B-layers. A- and B-layers are rotated by 60° around the c- 
axis with respect to each other. Flat and buckled Al-Si layers 
are indicated with and without an asterisk, respectively. 



Recent X-ray diffraction experiments on nH-C&A\Si 
single crystals revealed several stacking variations of two 
different types of Al-Si layers, denoted as A and B in the 
following (see Fig. [J). The A- and B-layers differ by a 
60° rotation around the c-axisi 13 ' 14 Besides the simple 
1H structure characterized by a stacking, two 

more stacking variants were found: 5H with a AABBB 
and 6H with a AAABBB sequence. Stacking of the A 
and B layers induces an internal stress on the structure, 
causing a buckling of boundary layers, i.e. layers with an 
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unlike neighboring layer. Subsequent investigations show 
that the superconducting properties strongly depend on 
the kind of stacking of A- and B-layers . 14 i 19 ' 20 ' 21 How- 
ever, all ab-initio calculations, and analyses of the ex- 
perimental data based thereon, so far have assumed that 
the stacking of Al-Si planes along the c axis is either uni- 
form, or completely disordered. To clarify the complex 
interplay of stacking variants, buckling, soft modes and 
superconductivity, further investigation of the e-ph prop- 
erties of different stacking variants of CaAISi are highly 
desirable. 

In this paper, we investigate the effects of pressure 
on the superconducting properties of nif-CaAISi by ex- 
periments on single crystals and ab-initio calculations of 
the e-ph properties. Previous experiments on polycrys- 
talline samples of CaAISi and SrAISi have shown that 
T c of CaAISi increases under pressure while T c of SrAISi 
decreases^ In an ab-initio study, Huang et al. have pro- 
posed that in liJ-CaAISi an enhancement of the e-ph 
coupling associated to the ultra-soft phonon modes can 
overcome the negative contribution of the other phonon 
branches and lead to an increase of T c under pressure. 
Such an ultra-soft phonon mode was not obtained for 
SrAlSi.i£ 

However, recent experimental and theoretical results 
question these findings. First, until now the pressure 
dependence of T c in polycrystallinc samples has been 
interpreted assuming a IH stacking variant. However, 
we will demonstrate here that the pressure dependence 
of T c of lH-CaAISi single crystals is markedly different 
from that of polycrystalline samples. Secondly, in these 
calculations^, the soft phonon branch in CaAISi is un- 
stable in some part of the brillouin zone (BZ), already 
at ambient pressure. This instability was not confirmed 
by recent ab-initio calculations and neutron scattering 
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Finally, Huang et al. assumed a uni- 



form compressibility, which is very unlikely to occur in 
layered materials. 

The aim of our analysis is two-fold: on the one hand 
we wish to investigate the interplay between soft modes 
and superconductivity in lif-CaAISi based on ab-initio 
calculations and experiments on single crystals. On the 
other hand, we wish to understand how stacking variants 
affect the pressure dependence of the superconducting 
properties of CaAISi. 

Single crystals of CaAISi were grown by the floating 
zone method. Details of the crystal growth and charac- 
terization are described elsewhere^ The superconduct- 
ing transition temperature (T c ) under pressure was de- 
termined by measuring the magnetic susceptibility in 
a SQUID magnetometer (Quantum Design). A Cu-Be 
piston-anvil-type pressure cell was used to apply quasi- 
hydrostatic pressures up to P ~ 10 kbar with silicon oil 
or Fluorinert as pressure transmitting medium. To mon- 
itor the pressure inside the pressure cell, we performed 
in-situ measurements of the T c of Sn (purity 99.999%) 
or Pb (purity 99.9999 

Figure [5] shows the temperature dependence of the 
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FIG. 2: (color online) Left: Temperature dependence of the 
susceptibility for CaAISi single crystals at different pressures. 
The numbers next to the data and in the parenthesis cor- 
respond to the applied pressure (kbar) and the sequential 
order of the measurement runs, respectively. The extrapo- 
lation method to determine T c is demonstrated with solid 
lines. Right: Pressure dependence of T c for (a) 1H-, (b) 5H-, 
(b) 6-ff-CaAlSi single crystals. For comparison we plot the 
previous data for polycrystaline sample from (Ref. fla ) in (c) 
and the theoretical curve for T C (P) in (a) for lff-CaAlSi with 
(red) dashed line.— 



magnetic susceptibility (x(T)) of single crystals of nH- 
CaAlSi under pressures up to ~ 10 kbar. At ambient 
pressure, the T c 's of IH, 5H, and 6-ff-CaAlSi are 6.50 K, 
5.95 K, and 7.89 K, respectively, consistent with previous 
reports^ The superconducting transitions are relatively 
sharp for IH and 6-ff-CaAlSi with a transition width of 
AT C ~ 0.1 K at 80% of the diamagnctic shielding, indicat- 
ing a good sample quality. For 5-H-CaAlSi, the transition 
is broader (AT C ~ 0.4 K), but still sharp enough to follow 
its pressure dependence. With increasing pressure, the 
superconducting transition for liJ-CaAlSi clearly shifts 
to lower temperatures without significant broadening of 
the transition. In contrast, T c in 5H and 6-ff-CaAlSi in- 
creases with pressure. For all samples, the T c 's and the 
shape of x(T) trace taken after gradually releasing again 
the pressure agreed with the data collected with increas- 
ing pressure. 

Our measurements reveal that the pressure variation 
of T c of n.ff-CaAlSi depends strongly on the stacking se- 
quence. For the detailed comparison, T c was identified 
as the temperature where the extrapolation of the steep- 
est slope of x(T) intersects the normal-state susceptibil- 
ity extrapolated to lower temperatures. Using different 
criteria to determine T c , e.g. the mid-point of the tran- 
sition, changes the absolute value of T c , but does not 
alter the relative variation of T c . Surprisingly, for IH- 
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CaAISi, T c decreases linearly with a rate of AT C /T C = 
-0.03 K/kbar, in contrast to a previous report on poly- 
crystalline samples.— For 5-H-CaAlSi the pressure depen- 
dence is slightly nonlinear with an initial slope of +0.013 
K/kbar, while 6i/-CaAlSi exhibits a more pronounced 
nonlinear behavior and saturation at T c ~ 7.95 K already 
at ~ 10 Kbar. 

In order to understand the effects of pressure and 
stacking variants on the superconducting properties 
of CaAISi, we performed ab-initio (DFPT) calcula- 
tion o 24 i 25 i 26 of the electronic and vibrational properties of 
n_ff-CaAlSi as a function of pressure. The computational 
details are the same as in Ref. [TH, unless otherwise spec- 
ified.— Besides the identified IH, 5H and 6H stacking 
variants, we also considered hypothetical 2H and 3H- 
CaAlSi, characterized by an |AB|AB| and AAB|,4AB| 
stacking of the Al-Si planes, respectively. 

In the investigated pressure range up to 100 kbar, we 
relaxed the internal parameters of all structures at inter- 
vals of 20 kbar. The details of the optimization, together 
with the detailed band structures and phonon disper- 
sion relations, will be published elsewhere. We subse- 
quently fitted the E vs. V curves to a Birch-Murnagham 
equation-of-state, to obtain the equilibrium volume Vq 
and the bulk modulus Bq for each system. All equilib- 
rium volumes and bulk moduli are close to each other 
(68.6 < V < 69.0 A 3 , 49.6 < B Q < 51.8 GPa). 
In sign and magnitude, our calculations correctly repro- 
duce the corrugation ("buckling") of the Al-Si planes in 
the 5H and 6H structures. A corrugation is also obtained 
for 3-ff-CaAlSi. In contrast to previous suggestions,— the 
corrugation in the Al-Si layers of the 5H and 6H phases 
is not reduced by applying pressure. Rather, it remains 
unchanged for 6H, or becomes even more pronounced for 
5 H- CaAISi. The compressibility is anisotropic, k c /k a = 
3.8, as expected for a layered material. At all pressures, 
the total energy of 2_ff-CaAlSi and 3iJ-CaAlSi is higher 
than in 1H, 5H and 6i/-CaAlSi, the latter being almost 
degenerate in energy.— 

Using the relaxed structures, we calculated by linear 
response the phonon dispersion relations w q .„ and the 
e-ph lincwidths 7 q .„ on a regular grid in q— space. The 
phonon dispersions and the phonon Densities of States 
(phDOS) were then obtained by Fourier interpolation of 
the dynamical matrices, while the Eliashberg spectral 
functions a 2 F{uj) and the total e-ph coupling A's were 
calculated by summing over individual lincwidths and 
frequencies. 

In Fig. [H we display the phDOS of niJ-CaAISi, calcu- 
lated at the theoretical equilibrium pressure (P = 0) and 
at the highest pressure considered (100 Kbar). At P = 0, 
the 1H, 5H and 6H stacking variants all show a peak at 
~ 7 meV, associated with the out-of-plane vibrations of 
the Al-Si planes. Under pressure, roughly one half of the 
corresponding phonon states soften and drive the system 
to a structural instability at P ~ 80 Kbar.— The hypo- 
thetical 2i?-CaAlSi and 3ff-CaAlSi, on the other hand, 
do not display the ultra-soft phonon peak at P = 0, and 
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FIG. 3: (color online) Linear Response Phonon DOS (phDOS) 
for different nH-CaAISi at P=0 and P=100 kbar; in the right 
sub-panel we show an enlargement of the low-energy region. 



all phonon modes harden with increasing pressure. 

At present, we cannot conclusively explain the different 
phonon softening behavior of the systems. However, by 
comparing our (partial) phDOS's obtained for different 
nH stacking variants with those of previous VCA cal- 
culations ji^iil we observe that phonon softening mainly 
involves Al out-of-plane vibrations. Phonon softening ap- 
pears if a block of three or more like planes is found in 
the sequence. This probably reflects the tendency of Al 
to form bonds with other Al atoms in the neighbouring 
layers rather than sp 2 bonds with Si in the same layer, if 
the c lattice constant is small enough. 

Our results strongly differ from those of Ref. [l6| for 
1H— CaAISi. At P = 0, our phonon frequencies are all 
real, and even at higher pressures only a few become 
imaginary, while in Ref. [l6| the whole out-of-plane branch 
becomes unstable. The most likely reason for this dis- 
crepancy is that we optimized the structural parameters, 
while the previous work assumed an isotropic compress- 
ibility. Also, a careful convergence of phonon frequencies 
with respect to k-sampling and basis set size is crucial 
to obtain the correct behaviour of the ultra-soft phonon 
modes. 
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FIG. 4: (color online) (a) Pressure dependence of the DOS at 
Ef and (b) frequency of the lowest-lying phonon branches at 
the H point for various nff-CaAlSi.— 
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If the coupling to the electrons is strong enough to 
overcome the opposing effect of the hardening of the other 
phonons and of the decrease of the electronic DOS at the 
Fermi level, a soft phonon mode can lead to an increase 
of T c with pressure. In Fig. [4] we show the calculated 
pressure dependences of the electronic DOS at the Fermi 
level, N(Ep), and of the lowest lying phonon frequency 
at the H point {uj h ) for IH, 5H and 6iJ-CaAlSi. We 
estimate the pressure evolution of the partial e-ph cou- 
pling constant associated to the soft phonon mode, using 
the Hopfield formula: A = N(E F )I 2 /ojfj, where / is the 
e-ph matrix element. Assuming that / is independent of 
pressure, we estimate that the partial A increases by ~ 
60% for 1H and 6H, and ~ 50% for 5H. For 2H and 3H- 
CaAlSi (not shown), N(Ef) decreases and ujh increases, 
resulting in a net decrease of A. From this argument, 
one could expect a very similar (increasing) behavior of 
T c with pressure for n = 1,5,6. This is, however, not 
observed experimentally. 

Figure [5] shows the evolution of the Eliashberg func- 
tion a 2 F(uj) of liJ-CaAISi for which we ran a full e-ph 
calculation as a function of pressure. In agreement with 
previous calculations , 12 ' 17 at P = there is a large peak 
at uj ~ 7 meV, associated to the Al-Si out-of-plane vibra- 
tions in the fc z =7r/c plane, implying that the correspond- 
ing phonon states have a considerable coupling to elec- 
trons. At higher pressures, following the behavior of the 
phDOS, this peak splits into two peaks, one of them soft- 
ening, the other hardening with increasing pressure. The 
pressure-induced transfer of the spectral weight from the 
hard to the soft part of the low-lying peak affects both the 
total e-ph coupling parameter A = 2 a 2 F(uj)oj~ 1 doj 
and the logarithmic-averaged phonon frequency (u>i n ). 



Using the Allen-Dynes equation^ 2 - 
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FIG. 5: (color online ) Evolution of the LDA Eliashberg func- 
tion (open circles) a 2 F(uj) and A(oj) with pressure (kbar), 
for Iff-CaAlSi. Red solid lines correspond to the simplified 
model with Gaussian peaks, described in the text. The re- 
sulting (uJin) and A are listed in Table. [I] The shaded area 
highlights the variation of energy and intensity for the ultra- 
soft phonon modes. 
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where fi* is the Coulomb pseudopotential, fixed to /Lt* = 
0.1 in the following (see Tab. [I]), we obtain a decrease of 
T c with pressure for lif-CaAlSi, in agreement with our 
experiments. A slight variation of /i* does not alter this 
general behavior. 

In stark contrast to the previous report for IH- 
CaAlSi,— the effect of softening of the ultra-soft phonon 
modes is apparently not sufficient to overcome the oppos- 
ing contributions from the other phonon branches. For 
further model calculations, we used a decomposition of 
a 2 F(ui) of liJ-CaAISi into Gaussians as shown in Fig. [5] 
Using this decomposition reveals that the behavior of T c 
with pressure follows that of the spectral weight of the 
ultra-soft phonon mode, which is decreasing with pres- 
sure in lff-CaAlSi. The change of the spectral weight of 
the ultra-soft phonon mode determines the behavior of 
T c under pressure, more importantly than the decrease 
of its phonon frequency. 

For 5if and 6ii-CaAlSi, due to their very large unit 
cells (15 and 18 atoms, respectively), and the extreme 
sensitivity of the results to computational parameters, we 
did not perform a full e-ph calculation. Although there 
are some subtle differences in the electron and phonon 
dispersions, the results shown in Fig. [3] and pimply that 
the pressure dependence is qualitatively very close to that 
of liJ-CaAISi, (i.e. the N(Ep)'s and ujh decrease with 
the same rate, and a structural instability happens at el- 
evated pressures). For this reason, we assume that their 
a 2 F(uj)'s have very similar characteristics to that of \H- 
CaAlSi. The different pressure dependence of T c for 5H 
and 6H therefore has to be attributed to a different (in- 
creasing) behavior of the spectral weight for the ultra- 
soft phonon mode under pressure. This behavior could 
reflect either small differences in e-ph matrix elements 
due to the buckling of some planes, or a different number 
of phonon modes that soften under pressure. 

Another aspect, which we cannot rule out completely 
is, that in 5H or 6-H-CaAlSi multi-band effects become 
relevant. Recently Lupi et al. reported an anisotropy of 
the optical response in the superconducting and the nor- 
mal states for a crystal with mixed 5H and 6H phases. 21 
The superconducting gap, determined by penetration 



P (Kbar) 


(win) (K) 


A 


Tc (K) 





139.2 (160) 


0.73 (0.60) 


5.35 (3.66) 


20 


139.9 (186) 


0.70 (0.47) 


4.86 (1.75) 


40 


134.0 (200) 


0.70 (0.43) 


4.65 (1.24) 



TABLE I: Calculated superconducting properties of 1H- 
CaAlSi as a function of pressure. T c was obtained by Allen- 
Dynes formula, with fi* — 0.1. For comparison, we give in 
parentheses the results for 2H. 
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depth measurements, also shows a sizable anisotropy, de- 
pending on the stacking sequences^ Further studies are 
needed to clarify this point. In particular, understanding 
the effect of buckling of the Al-Si layers on the anisotropy 
of the e-ph coupling appears to be the crucial issue. 

We would like to emphasize once more that the hypo- 
thetical systems 2H and 3ff-CaAlSi show a completely 
different behavior. First, they are energetically disfa- 
vored with respect to existing stacking variants, because 
of a non-optimized energy balance between the formation 
of AB interfaces and buckling. Second, they do not dis- 
play any soft phonon modes, which only appear if three 
or more Al atoms arrange in a sequence along the c-axis. 

All phonon modes harden with pressure, leading to 
an increase in (win), a fast decrease of A, and a net de- 
crease of T c with pressure with a rate of dT c /dP ~ -0.05 
K/Kbar. For comparison, the results for 2_ff-CaAlSi are 
also listed in Table HI 2H and 3i/-CaAlSi, therefore, do 
not represent proper models for other stacking variants. 

In summary, we demonstrate that for single crystals of 
n_ff-CaAlSi the behavior of T c under pressure depends 
crucially on the particular stacking sequence and the 
buckling of the Al-Si layers. T c decreases for 1H, and 



increases with pressure in 5H and 6-ff-CaAlSi. Previous 
experiments on polycrystaline samples gave results that 
are closer to the behavior of 6f/-CaAlSi, but adverse to 
the behavior of liJ-CaAlSi. One may speculate that the 
polycrystalline sample consisted mainly of phases with 
stacking variants other than li?-CaAlSi. Based on our 
ab-initio calculations we find a gradual softening of an 
out-of-planc phonon mode under pressure for n_ff-CaAlSi 
(n =1, 5, 6), which leads to a structural instability at 
higher pressures. In 1H— CaAISi, the softening is not 
strong enough to lead to an increase in T c , in contrast to 
previous calculations, while it is likely that this softening 
leads to the increase of T c with pressure in single crystals 
of 5H and 6#-CaAlSi. 
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